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For some purposes, particularly in connection with the study of the random structure 
of the lower ionosphere, using very low frequencies, it is necessary to find the detailed 
statistical properties of a random signal diffracting in free space. Mathematical tools for 
evaluating these parameters have been developed, and are applied in this paper. Allowance 
is made for the effect of sphericity of the wave incident on the ionosphere, and anisotropy 
of the irregular variations of signal is permitted. The case of oblique incidence of a wave 
on the ionosphere is also considered. 

1. Introduction 

There has of Lite been an increasing interest in the irreguhir structure of the lower iono- 
sphere, particularly in connection with scatter propagation. In order to interpret irregular 
fluctuations in radio signals observed at the ground in terms of irregularities in the ionosphere, 
it is first necessary to find how much the nature of the signal changes as it propagates from the 
ionosphere to the earth. This propagation, which takes place in free space, will be termed 
^^diffraction," in tlie sense used in a previous paper [Bowhill, 1957]. 

Most of the methods so far proposed for dealing with this statistical problem assume that 
the signal at the layer is like the signal produced when an infinite plane wave passes through 
a thin diffracting screen. The field at the screen is tlien resolved into a number of plane waves 
propagating in various directions, their amplitudes being given by the Fourier spectrum of 
the field variation. The signal at any distance from the screen is then given by the vector 
sum of the signals from tlie various plane waves, combined witli the appropriate phases. 
Booker, Ratcliffe, and Sliinn [1950] stated this principle, but Hewish [1952] was the first to show 
quantitatively how the magnitude of the irregular field variations changed with distance from 
the screen. He assumed a rectangular distribution of spatial wave numbers in the screen. 
Feinstein [1954] has used a rather dift'erent approach based on Huy gens' principle, which has 
also been exploited by Chernov [I960]. Ratcliffe [1956] has summarized the physical prin- 
ciples involved in these approaches. 

Pitteway [1958] has recently given an analytic method for calculating the field due to any 
particular assembly of irregularities; this, however, is difficult to apply in a statistically random 
medium. 

It is important to distinguish at this point between the various t^^pes of random screen 
postulated. The two types most commonly considered are screens which modulate the phase 
or the amplitude of the signal. A phase screen is said to be ''deep'' or ''shallow'' depending 
on whether the phase excursion of the signal is substantially larger or smaller than one radian. 
The deep phase screen problem has been approached by Hewish [1951] by using a sinusoidal 
variation of phase with distance, but the extension of a detailed theory to this case presents 
considerable diflftculties. The shallow phase screen, on the other hand, is exactly equivalent 
to a shallow amplitude screen [Bowhill, 1959]. This interchangeability of phase and amplitude 
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is exploited in the phase-contrast microscope, in which variations in pliase are changed into 
variations in amplitude by artificially retarding the phase of the specular component by x/2. 
It is well know^n that this only gives an improvement in visibility when the pliase variations 
are small (i.e., the phase screen is shallow). 

The results of the analysis in this paper will be applied principally to the random diffrac- 
tion of very low frequency radio waves. The fading for these frequencies is usually shallow 
indicating that a shallowly modulated random field at the ionosphere must be responsible. 
Attention will therefore be confined to shallowly modulated screens. 

To illustrate the method of analysis, it will first be applied, in section 2, to a one-dimen- 
sional random screen. The extension to two dimensions is made in section 3. Section 4 deals 
with the correlation between amplitude and phase in the signal, and section 5 with the effect 
of wave sphericity. The work is extended to oblique incidence in section 6. 

2. One-Dimensional Diffraction Problem 

In treating the spatial correlation function and spatial frequency spectrum of the signal 
variations, it is much more convenient to use purely real quantities, rather than the complex 
variables of amplitude and phase. This enables the use of complex Fourier transforms to re- 
late field variations over the screen to the complex angular spectrum of plane waves which 
they produce. 

For an amplitude screen, which alone will be considered (in view of the equivalence 
indicated above), the electric field E at the screen is given b}^ 

E{x,,0)=A{x,) 

where A(:Xq) is a purely real quantity. The field at a distance z from the screen is then defined 
to be 

E{x, z) = [Ei{x) + iE2{x)] exp (27ris/X) (1) 

where Ei{x) and E2{x) are both purely real, and represent the phase components of E(x,2) in 
phase and in quadrature with the undiffracted portion of the incident plane wave. For a 
shallow screen, Ei(x) gives the amplitude variation, and E2{x) the phase variation. In the 
analysis which follows, various statistical parameters of the quantities Ei(x) and £2(^0, which 
are actually observed in experiments at ground level, are evaluated in terms of the corre- 
sponding parameters of the original amplitude variation A(xo). 

Using the notation of Booker, Ratclifte, and Shinn [1950], we have for the angular 
spectrum at the screen, in terms of s = sin 6, where 6 is the scattering angle, 

J 00 
A{xo) ' exp {—2'7risXo/\) dxo (2) 

where X is the wavelength. The diffracted field E(x,z) is given by the Fourier transform of 
this, together with the pliase factor exp (2Tricz/\), where c = cos d. 

E(x,z)=- exp {2Tri(sx+cz)l\}ds * A(xo) - exp {—27risXo/\)dxo 

'^ t/ — 00 %J — m 

=- ^(xo) • exp {2'Ki[{x—x^s+cz]l\]dXi^ds. (3) 

By the definition of Ei{x) in equation (1), 
£'i(x)=:;^/^ j j A{x^) . exp {2Tri[{x-XQ)s-{\-c)z\l\]dx4s\ 

=2x I j ^(^0) [^^"^P [2Tri[{x-x^)s-{\-c)z]l\] + exp{-2Tri[{x-x^)s-(\-c)z]l\}]dx4s 
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siiicc^ ^ I (./•()) is piirel}^ real. The spatial spectrum Pi{v) of Ei{x) is now defined by a similar 

(wprc^ssion to (2) : 

/^W=J " E,{x) . exp (-27ri.^:/X)rf^-^ ^^^^ A{x,) • [exp [2^i[^x-x,)s-vx- z{l-c)]l\] 

+ exp { — 2'Ki [ {x— Xo) s + vx—2(l — c) ]/X } ]dxodxds. 

Since the integrand is a well-behaved function, and all tlu^ limits of the integration are 
— oo to + ooj the integration can be carried out in any order. In fact, it is convenient to in- 
tegrate with respect to first s and then x. With the approximation (vahd for small s) that c = 

1—8^/2, Pi{v) is easily evaluated as 

PiM=M ^ ^(a-e)[oxp {27r/(-.Jo+,e.72)/X}+exp {2iri{-vx,-zv'l2)l\]]dx, 

= COs(7r2^7^) • A{'^()) ' t'xp (—2TrivXQ/\)dXo 

*/ — 00 

or 

AW=AW •cos(xj.7X) (4) 

from e(iiiali()ii (2). Tliis ecjuatiou gives th(* spatial sjXM'lrimi of E\(x) at any dislancc^ z from 
the screen. TJie power spectrum TT^(j^) of this variation is defincMl by 



W,(v)=P,{v) ./1=W (5) 

the asterisk denoting the complcw conjugate, and the bar denoting averaging of tlie result 
'^over systems'' — i.e., tlie mean for a mnnber of statistically similar screens. This quantity 
is essentially real, and is a relatively smoothly varying function of v, because of the process of 
averaging. Tt can be related by (Hpiation (4) to the spatial powcM- sjXM'trum \\\)(i') at tlu^ scrcHMi : 



W,{V)=P,{V) . Pt{v) • cos' [TTZV'/X) 

or 

W^(v)=Wo(v) . cos' (tzv'/\) (6) 

Avhere 



Similar quantities P2(p) and W^iv) can be defined for tlie quadrature component E2 of the field 
(see equation (5)); a similar analysis to the preceding gives 

P2(v)=Po{v)'Sm(T2vy\) ^ 
H^2 W =Wo{p) •sin2(7r.3^7X) / ^^^ 

Equations (6) and (7) enable the spatial correlograms of Ei(x) and E2(x) in the diffrac- 
tion pattern to be determined. Assuming that the correlogram of A(x) is given, the Wiener- 
Khintchine theorem can be used to give Wo{p). The power spectrum Wi(v) is then calcu- 
lated from equation (20), and the coirelogram of Ei(x) found by the inverse Wiener-Khintchine 
theorem. 

This one-dimensional analysis, in which the random variations in signal are taken to 
occur in the 2:-direction only, represents an unrealistic physical model. However, as some 
existing work [Hewish, 1951, 1952; Jones, Millman, and Nertney, 1953] has assumed this 
type of variation, some results worked out for this case will be quoted in the next section, 
for comparison with the two-dimensional results. 
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3. Two-Dimensional Diffraction Problem 

In this section the diflfraction pattern of a two-dimensional random screen is considered 
using an analysis closely similar to that employed in the previous section for the one-dimen- 
sional case. The modulation is once again assumed to be in amplitude only, as there is again 
an equivalence between shallow phase and amplitude screens. The effect of anisotropy of 
the screen will be taken into consideration by defining a spatial power spectrum Wo{v,x)j 
related to a two-dimensional spatial correlogram PQ{a,P) of the field A(x,y) at the screen by 
a two-dimensional form of the Wiener-Khint chine thereom. As the field A(xo,yo) is purely 
real, the correlogram po{a,P) is given by 



Aix,y) A(x+a,y+p)-A(x,yy 
A'ix,y)-A{x,y)' 



The field A{x,y), the mean field over the whole screen, has a constant phase at all points since 
Aix^y) is real, and propagates as a single plane wave. It is convenient to imagine it removed; 
A(x,y) then has both positive and negative values, and will be defined to have a mean square 
value of unity; so 




. (8) 

A\xo,yo)-- 

and the correlogram is given by 



Po{a,P)=A(xo, yo)A(xo+a, y^+P). 

The single plane wave, which will be called the ^'specular component" (by analogy with 
specular reflection from a mirror, compared with irregular reflection from a rough surface), 
can be reintroduced when necessary into the diffraction pattern as a constant added to Eiix^y). 

Proceeding exactly as in section 2, the two phase components of the diffracted field are 
defined by the vector field at (x,y,z) : 

E(x,y,z) = [E,(x,y)+iE2ix,y)] exp (27riz/\). (9) 

The phases of the diffracted waves at z are given by purely geometrical considerations. Let 
the plane of one wave intersect the {x,z) plane at an angle sin" ^5, and the {y,z) plane at an 
angle sin~V. Then the angle between the wave front and the (x,y) plane is 

tan-i{ (^-2_ i)-i_^ (,^-2_ i)-i }i/2 
and the phase of the wave at a distance z from the screen is 

If s and r are small, this can be expanded to give, to a first order, 

2iriz{s^/2 + ry2)l\ 
The corresponding equation to (3) for E(x,z) is therefore 

E(x,y,z) exp (—2Triz/\)y 

=- I I I I A(xo,yo)'exp [27ri[{x-Xo)s+(y—yo)r-s^zl2—r^z/2]l\}dsdrdxodyo (10) 
and by similar steps to the previous analysis, it can be shown that 

Pi{v,X)=Po(v,X)'COs{irz(p' + X')l\} (11) 
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Avlit're /•'u(v,x) is the two-diiiionsional spectrum of tlic field at the screen. Similarly 

and the spatial power spectra of Ei(x,y) and E2{x,y) are given by 

W,ip,x)=Woiv,x).cosnTz(v'+X')/\] '] 

Y (12) 

These spectra are proportional to the area of the diffracting screen; this will be shown below 
in deriving (16). The spatial correlogram pi(a,/3) of Ei{x,y) is given by tlie usual relation 



E,{x^y)E^^+a^ 
E\{x,y) 

:since the mean value of Ei{x,y) is zero. On comparing this definition with (9), it will be 
.seen that (13) has the term Elix,y) in the denominator, whereas, from (8), A\xo,yo) = l. In 
fact, the mean square fluctuation in A{xo,yo) at the screen appears, at some distance away, as 
a fluctuation partly in E2(x,y). It is easily shown from (12) that 



EKx,y) +EI{x,y) =A\xo,yo) = 1, (14) 

the sum of the mean square values of Ei(x,y) and E2{x,y) r(Miiainini2." constant at all distances 
from the screen. These mean square values represent, in the case of an amplitude screen, the 
amount of fluctuation of the amplitude and the pliase respectively at a distance z from the 
screen. 

The autoco variance of Ei can be calculated from the spectrum Wi(v,x) in (12) by the 
Wiener-Khintchine theorem; E\{x,y) is then its value for a=^=Q. Before this can be cal- 
culated, the power spectrum Wi{v,x) of the signal A{xQ,yQ) at the screen must be found. Some 
assumption is needed for tlie spatial characteristics of thc^ signal. They will be assumed to 
be as follows: 

i. A{x^y) has a Gaussian correlogram in all directions. 

ii. The contours of the two-dimensional correlogram are elliptical. 

iii. The ellipses have their principal axes oriented along the x and y axes of the co-ordinate 
:system . 

The p=0.61 contour on the two-dimensional correlogram is used to define two structure 
sizes Di and 7>2 of the screen in the x and y directions. The correlogram of A{x,y) is then 

Po(a,/3) =exp{ - (a72P?) - (^72Z>i) } . (15) 

The angular spectrum at the screen is given by an equation like (2) 

Po(s,r)= J I A{xQ,yo)-exp { — 2Tri(sxo+r7jQ)/\}dxodyo 

and the power spectrum is given b}^ 

Wo{s,r) =Po{s,r) 'Pt(s,r) =mean of ^ I I A(xo,yo) • exp [—27ri(sxo+ryo)/\}dxodyo 

X j I A{xo,yo) • exp { 2Tri(sxo+ryo)/\}dxodyo V 

A(xo,yo)-A{xo,yQ)'exp{2Tri[(xo—Xo)s+ (yo—yo)r]/\}dxodyGd4dyo 
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=m: 



where Xq and y'o are variables of integration, analogous to ^^o and yo, but differentiated from 
them to avoid confusion. So 

Wo(s,r)= I I I A{xo,yo)'A(xQ+a,yo+P).exp [2Ti(as + MMd^odyo(Ia djS 
where a=Xo — Xo, l3 = yo — yo. By (14), it is evident that 



A{xo,yo)-A{xo + a, 7/o+0) = Po(o;,/3) 
and therefore 

Wo(^7)= JJJJ Po(a,P)-exp{2TiiaS + MNdxodyodadp, 

Now neither Xq nor y^ appears in the integrand; so the integration with respect to Xq and y^ will 
give a Wo(s,r) which is proportional to the total area >§ of the diffracting screen: 

Wo{s,r) = S- \\ po(a,l3)-exp[2Ti(as + pr)l\}dadl3. (16) 

Substituting for po{a,l3) from (15) and integrating, 

W'^o(5/) = 27rSAA.exp{-27r^(AV+AV)/X}. (17) 

The spatial power spectra of Ei{x,y) and E2{x,y) are found by substituting from (17) into (12): 

T^i(.,x)-2xSAA.exp{-2^XAV+AV)/V}.eos2{^^(.2+X^)/X}\ , 

TT'2(^,x)=27rSAZ^lexp{-2^^(AV + AV)/X^}.sin^{^^(^^ + X^)/X}/ ^ ""^ 

It is now possible to find the numerator of (13); the covariance oiEi(x,y) and £'i (x + a^y^^). 
The inverse transform corresponding to (16) is 

E,{x, y) E,{x+a, ^+/5):=^, JJ W,{v, X) • exp{-2Ti{va+X^)/\}dv dx (19) 

with the value of Wi(v,x) given by (18). The integral is evaluated by putting 

cos^T2{p'+x')l\}==l/2+(ll4)exp{2Triz(p'+x')l\] + (ll4:)exp{-2iriz(p'+x')l\}. (20) 

As (19) is the inverse of (16), the term 1/2 gives simply (l/2)po(a,l3). The two remaining 
terms in (40) are identical except that the sign of z is reversed. The integral can be simplified 
further by separating completely the integrals with respect to v and X. This gives 

E^{x,y).E,{x+a,y+l3)^'^^ {I{a,D,,z)I{&, D„ z)+I(a,D„ -z)I{p,D,, -s)} + (l/2)po(«+/3) 

(21) 
where 

I{a, D,,z)={ exp { 2Tvi[av+ {z +xiZ>?/X) v']/\ } dv 

*J — CO 

_X(l+a|^ exp {-(1/2) tan"' ai-{l-ia,)a'l2D\{\+a\)] (22) 
V27rZ/i 

where the important substitutions 

\z \z 

have been made. These parameters will appear as the measures of distance on all the graphs 
depicting the behavior of the various statistical parameters of the diffraction pattern. 

280 



Finnlly, on substituting from (22) aud (21), 



/':,{x,y)-h\{x + a,y+(i) = (l/2)exp{-ay2D{'~l5y2D,'} 

+ (l/2)(l+a,^)-"^(l+a/)-i'^cxp{-aV2A-(l+«i^)-^/2Z)i(l+aD} 
Xcos{aia72A'(l+ai') + a2(3V2A.'(l+a/(1+ffl22)-(t!in-V/i + t!in-'«2)/2}. (23) 

By putting a = /3=0 



^1^(^,2/) = 1/2+ (l/2)(l + (i,^)-"*-(l+a/)-"*-cos{(tan-'a, + taii-'a2)/2} 



(24) 



The correlogram pi{a,P) is given by the ratio of (23) to (24). Tlu^ expressions lor the auto- 
covariance of E2(x,y) are identical, except that the positive signs, (UMioted ahovc^ with asterisks, 
become negative. Some special cases will now be investigated, 
(i) Near the screen ai = a2 = 0, and 



E,{x,y)'E:{x+a,y+^) = oxp(-a'l2D{'-^'l27l/) = Po{a,^) 



E^(j',y) = i,K2'^;y) = 0. 

as w^ould be expected, since Ei(x,y) = A^Xo^yo) at the screen. 

(ii) Very far from the screen ai ,2-^ °° , tan ~ ^ai = tan ~ ^ao = t/2 , 

and 

E,{x,y)E,{x + a,y + i:i) = (l/2)oxp{-a'/2D,'-l3y2l)2'} 



E,\x,y) = K,'(:x,y) = \l2 
p,{a,p) = p,{a,p) = cxp[-a'i2/K'-^'l2D2'} 



(25) 



The second result in (25) means that the space correhition of both amphtuch^ and phase at a 
great distance from an amplitude screen is just the same as the correlation at the screen. This 
is not true at intermediate distances, as we shall shortly see. The first result in (25) sliows that 
the fluctuations occur equally in amplitude and phase unchM* these conditions. The phase- 
amplitude locus of the instantaneous signal at largi^ distances th(M'(*r()r(^ forms a cii'cular ''Cloud'' 
around the vector denoting the mean carrier signal, as in figure 1. 
(iii) Isotropic screen Dx = T)y=E>^ and ai = a2 = a. 



E,{x,y)E,{x + a,y + p) = (l/2)cxp{-(a'+^')l2n'} 

+ (l/2)(l + a2)-i/^.exp{-(a2+/32)/2D2(i + a2)}.cos{a(c.2+/3^)/2/^2(i+a2)_tan-ia}. (26) 

The quantities a and /3 appear in this expression only as (a^+iS^), so the contours of correlation 
are circles at all distances, and the diffraction pattern is never anisotropic. Also 



^?(x,y)-l/2+(l/2)(l + a^)-i 



LOCUS OF RANDOM 
COMPONENT AS ONE 
AXIS IS TRAVERSED 




SPECULARLY 
REFLECTED 
SIGNAL 



Figure 1. A possible phase — amplitude diagram for 
the diffracted signal a long way from a random 
screen. 



^ E,(x,y) 
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the asterisk signifying, as previousl}^, the sign which is to be reversed to give the corresponding 
value for Ef(x~y). Ei^{x,y) and lE^(x~y) are plotted on figure 2. It can be seen that the mean 
square values of Ei and E2 are virtually equal at distances from the screen greater than that 
corresponding to about a =4. It will be shown later that, under these circumstances, the 
diffraction pattern is not correlated with the signal at the screen. 
Slight rearrangement of (26) gives 



E,{x, y) E,{x+a, 2/+/3) = (l/2). exi^{-(a'+P')l2D'} 

+ (ll2)(l + a')-Kex^{-{a'+p')l2D\l + a')} [cos{a(a'+p')l2D\l+a')} 

+ a sin {a{a'+^')l2D'{l + a')}l (27) 



(iv) One-dimensional screen 
From (24), 



D, = D; A-^«>. Soa2^0. 



Ei\x, :?/) = l/2 + (l/2)(l + a2)-i/4.cos{(tan-i a)/2} 
= ll2 + {ll2){[l + {l+arm+aT'l2} 



1/2 



This also is plotted in figure 2, for comparison with the curve for the two-dimensional case. 
Evidently the mean square values of Ei{x,y) and E^ix^y) become nearly equal much closer to a 
two -dimension ally irregular screen than a screen which is irregular in one dimension only. 
This illustrates that a one-dimensional model may give a quite inadequate representation of two- 
dimensional diffraction phenomena. For this case, 



E,{x, y) E,{x+a, y+p) = {ll2) exp(-aV20^) 

+ (l/2)(l+a2)-i'4exp 
— (tan-ia)/2}. 



-a'l2D'\l+a')}'COs[aay2D'{l+a') 




Figure 2. Variation of mean square diffracted signal 
with distance from two kinds of random screen. 
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It will !)(' noted from (27) that the correlogram is not Gaussian in form, though Ei{x,y) has a 
Gaussian correlogram for a=0, and both Ei(x,y) and E2ix,y) have Gaussian correlograms as a 
-^00 (see equation (25)). This introduces some difficulty as to what is to be called the ^'struc- 
ture size" d of the diffraction pattern, for comparison with experiment. The following defini- 
tion is adopted 

c/ = Structure size in x direction = ( — d-p/^c>:^Oa=^//-() 

= Di for the correlogram po(a:,i3). 

This definition bases the structure size on the Gaussian cnrve which is the best fit to the correlo- 
gram, near the origin of a. It can best be found by expanding (23) in powers of a^ and /S^, and 
neglecting terms of order higher than the second; using equations (13) and (24), 

l + (l+aD"'^'(l+ai)-i/4cos {(tan-^ai + tan-^a2)/2} 
'' '^l + (l +«?)""*(] +al)'"' [cos (tan-' ai + tan-' a2)/2-ai sin {(tan"' a,+tan-' a2)/2]] 

l + {l+af)-"*{l+al)-"* cos {(tan-'a, + tan-'a2)/2} 
'/ ' ~ l + (l+a?)-"''(l +al)'"' |cos (tail-' a. + tan"' a,)/2-a, sin {(tan"' a, + tan-' a2)/2}] ' 

(28) 

To help in imdcrstandiiig the significance of these results, a number of special cases are 
examined below. 

(v) Near the screen, Ut and a2-^0, and di-^Di, (h^-lh for Ei{x,y), as found previously unch-r (i). 
For the quadrature component E2ix,y) however, it is easily shown tliat 



f K'+2Krs+l Y" 
'~^ \ 5K'+2K'+l f 



(29) 

d,. 



where K^D^jD.. 

(vi) Isotropic screen Di = D2=D, K=l. 

Near the screen, d=D for Ei(x,y) as found previously. 

From (29) d=DI^ for E2{x,y). 

This result implies that the structure size of the shallow phase variations near an amplitude 

screen (or of the amplitude variations near the phase screen) is 1/V3 times the structure of 

the screen itself. In addition, the correlogram of E2(x,y) is far from Gaussian; it can easily 

be shown by expanding (23) in powers of a that, near the screen, 

P2(a,p) = exp{-{a'+p')/2D'}.{l-(a'+p')ID'+{a'+^'mD'} (30) 

exactl3\ This is plotted in figure 3, and has an oscillating form. At any distance from an 
isotropic screen, (28) reduces to for Ei(x,y): 

''~ \ a*+a'+2 / 
for E2{x,y) : 

These relations are shown plotted in figure 4. The structure size of Ei{x,y) actually has a 
maximum at intermediate values of z. 
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Figure 3. Correlogram of the shallow phase variations 
near a two-dimensional amplitude screen. 



Figure 4. Variation of structure size d with dis- 
tance from a two-dimensional amplitude screen. 
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(vii) One-dimensional screen Bi^D) D2-^^, a2-^0, K->0. Near the screen, 

d=D for Ei{x,y), d=Dl^d for E2(x,y), 

showing the effect described in (vii) to an even greater extent. At any distance, 
for Ei{x,y): 



for E2(x,y): 






l + (l+a')"^^'cos{(tan-^a)/2} 



+ {\+a')-'/' lco^{{tan-'a)/2}-a sm{(tSin-'a)/2}] 



} 



1/2 



l-(1+a^)-^/^cos {(tan-^a)/2} \'/' 

-(l+(^')-'/Mcos{(tan-^a)/2}-asin{(tan-^a)/2}] / 

In the general case of an anisotropic screen, (28) can be used to find the structure sizes di 
and 6?2 in the x and y directions. To clarify the kind of behavior that may occur, figure 5 has 
been prepared. This shows how the contours of the correlation ellipse change shape with 
increasing distance from a random amplitude screen. The two screens shown are 

(a) isotropic, with A = 7^2=1 

(b) anisotropic, with Di = 2, Z>2=1. 

The distances from the screen were chosen as 0, oo , and TrZ^f/X. It is interesting to note that 
the phase correlation ellipse for the anisotropic screen does not keep the same shape as the 
screen is approached — its minor axis is reduced by a greater factor than its major axis. This 
is because as it approaches the limiting case of an anisotropic screen — namely, a one-dimen- 
sional screen — the major axis of the phase correlation ellipse near the screen approaches the 
true structure size Di of the screen, while the minor axis is reduced by the factor V5 established 
above in (vii) . 

The reason for dealing with this topic in such detail is to demonstrate that the structure 
size of the irregularities in the diffraction pattern produced by a shallow phase or amplitude 
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E|()c,y) 



ANISOTROPIC 
SCREEN, D| = 2, 



PHASE 
E2(x,y) 











Figure 5. Contours of the spatial correlogram of the 
wave diffracted from tivo kinds of amplitude screen, 
at different distances. 



inodulated screen is not always the same as that of (he sei'een itseH'. 'V\\c preceding lengthy 
mathematical argument is necessary to establish ihe pi'ecise effect of the diffraction process, so 
that existing (^xjXM'iiniMital (h-ita may h(^ int(M-|)reted with conlidence. 



4. Cross-Correlation Relationships for £i and E2 

A good deal of attention has been paid by various workers (lor instance, Jones, Millman, 
and Nertney [1958]) to the question of whether the signal in the (Uffraction pattern is the same 
as the signal at a random screen directly opposite to it. It must be emphasized that the 
question is not wliether the signals are statistically similar, i.e., have the same structure size or 
fading time, but whether the detailed changes in signal are similar at the two planes. Expressed 
in more exact language, the cross-correlation is to be evaluated between, for instance, the 
amplitude of the signal at the screen and the amplitude at a distant plane. 

Another cross-correlation coefficient of considerable importance is that between the 
amplitude and the phase of the diffracted signal [Bowhill, 1957]. Jones, Millman, and Nertney 
[1953] assumed tliat tliis correlation arises as a result only of correlation between the amplitude 
and phase of the signal at the screen. It is shown here, however, that they may be correlated 
even for a screen which changes the phase only or the amplitude only. 

The theory is developed for a shallow amplitude modulating screen, and is extended to the 
case of a phase screen. The quantity Ei{x,y) therefore represents the amplitude of the signal in 
the diffraction pattern (with the specular component removed) and E2{x,y) represents the phase 
of the signal. The quantities 



A{x,y)E^(x,y) and E2{x,y)E^{x,y) 



(31) 



are to be calculated; the known values of E\{x,y), El{x,y) and A\x,y) then give the cross- 
correlation coefficients pn between A(x,y) and Ei{x,y), and pi2 between Ei(x,y) and E2ix,y). 
The quantities (31) are computed from the relations 



^(^,2/)^i(2:,2/)=(l/2)^f(^,^) + (l/2)/LH.r,7/)-(l/2)[^(x,^)-£;i(^,7/)]^ 



^A(^,7y)^2(^,,?/) = (l/2)/^l(^,I/) + (l/2)AX^,^y)-(l/2)[^i(x,^)-^2(^,y)]^ 



■(32) 
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The last terms in these equations represent the areas under the power spectra of 
[A(x,y)—Ei{x,y)] and [Ei{x,y) — E2{x,y)], since the mean square value of any random variable 
with zero mean is given by the integral of its power spectrum with respect to frequency. 
The power spectrum of [A(x,y)—Ei{x,y)] can be found from (11); 



P,{v,X)^Poiv,X) . cos {7rz(v' + X')l\}. 
The spectrum of [A(x,y)—Ei{x,y)] is therefore 

Po{v,X)-P,(v,X)=Po(i^,X)[l-COs{irz(v' + X')l\]] 
and its power spectrum is given, exactly as in (6), by 

Wo(v,X)[l-COS {7rz(v' + X')l\}Y = 4 [W2{v,X)U-[W2{v,X)], 
and substituting in (32), 

Z(^;^)i\(^^(i/2)[:^f(^ 



(11) 



2/2- 



a relation which holds for any autocorrelation function at the screen. Using El{x,y) from (24) ^ 
the correlation pn is given by 

_a/2 • (l+a^/4)-^/^(l+a^/4)-^/^ cos {[tan"! (ai/2) + tan-^ ((^2/2)]/2} 



Pn=A{x,y)Ei{x,y)IE\{x,y)-- 



{l + (l+a?)-i/^(l+aD-'/'cos{[(tan-Ui+tan-ia2)/2]j 



for a screen with a Gaussian correlogram. 

For an isotropic screen, a^ = a2 = a, and this reduces to 

Pu={(l+a^)/(l + a72)}^'^(l+a74)-'. 

This function is plotted as the full curve in figure 6. The correlation is perfect near the screen^ 
since then Ei(x,y)=A{Xjy). It decreases progressively to zero with increasing distance^ 
reaching a value of 0.1 when a=7.2. 

The corresponding function for a one-dimensional screen, ai = a, ^2=0, is given by 



Pii=l 



2(l+aV4)-^/^ cos {[tan-^ (a/2)]/2} 
l + (l + a2)-i/^cos{(tan-iai + tan-ia2)}}^/ 



This function is plotted as the broken curve on figure 6. The correlation decreases much less- 
rapidly with increasing distance from the screen than for the two-dimensional case. 
To find the correlation pi2, (32) is rearranged to give 



E,ix,y)E,{x,y) = m-{l/2)[E,{x,y)~E,(x,y)Y 



ONE- DIMENSIONAL SCREEN 



TWO-DIMENSIONAL SCREEN 



Figure 6. Variation of the correlation of the diffrac- 
tion ^pattern amplitude with the amplitude at the 
screen as a function of distance for one- and two- 
dimensional screens. 
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and Uu; power spectrum of [Ei{x,y)—E2(x,y)] is given froni (11) bj^ 



[P: {v,X) -P, {v,X) ] [Pf (.,X) -I't {v,X)] 



So 



E,{x,y)E2{x,y) = {l/2\'S) ^f Wo{p,x) sin {27rz{v^+x^)/\}dvdx 

= (7rAA/2iX2) [[ [exp {2Tizii''+x^)/\}-exp{-2Triz(v'+x')/\}] 

•exp {-2ir2(Z??;'24.2}|x2)/X2}(;^(^X. 
= (7rAA/2iX2) J exp { 2Triv\z+TriD\l\)l\]dv exp { 27rix''{z+TriDl/\)/\]dx 

- r exp {27^^J/2(-^+7^^Z>?/X)/X}r/^ f exp {27^^x2(-^+7^^Di/X)/X}d:x\ 
and on inlegraiiiig this expression by the usual methods 



Also, from (24) 



4 ' Ei{x,y) ' I^{x,y) = l- (l+al)-'^\l+a^-'/' cos^iUm-' a, + t^n-^ a.;)/^^ 

and the correlation pi2 between Ei(x,y) and E2(x,y) is given b}^ 

_ (l+ai)-^/'(l+ai)-^^'sin{(tan-^ai + tan-2a2)/2} 
^''~[l-(l+a?)"^/Hl+tt2)~'^'cos2{(tan-iai+tan-ia2)/2}]i/2* 

For an isotropic screen, ai=a2=a; this reduces to 



(33) 



Pl2=(l+«^)-^/^ 

which is plotted in figure 7. The correlation between Ei{x,y) and E2{x,y) has a maximum 
value of 1/V2 at the screen, and decreases to zero with increasing distance. For an aniso- 
tropic screen, (33) must be expanded in powers of ai and (72. The correlation as ai and 02— >0 
is given by 

ai+a2 _ l+K^ . . 



Pl2- 



[3af+2aia2+3ai]^/2 [^K^J^2K^+^]' 



where K=DilD2 expresses the degree of anisotropy of the screen, as in (29). Whatever the 
value of K, (34) shows that this maximum value of the correlation pi2 lies between the values 

Pj2=1/V2 for an isotropic screen ai=a2 




Figure 7. Variation of amplitude-phase cross-cor- 
relation with distance from a two-dimensional 
amplitude screen. 
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and 

Pi2=l/V3 for a one-dimensional screen ai:$>a2. (35) 

The cross-correlation coefficients pn and pu between the phase components Ei{Xyy) and 
E2(x,y) must now be related to the amplitude and phase characteristics of the diffraction 
pattern. For an amplitude screen, it is apparent from figure 1 that for a shallow screen, 

Amplitude in pa ttern=^ (specular signal) +£'i(x,?/) 
Phase in pattern =^—£'2 (a:,?/) /(specular signal) 

since the positive phase direction is opposite to the positive E2 direction, taking an increase 
in phase to mean a retarding of the wave. This means that: 

Correlation between amplitude in pattern and at screen = pii, and correlation between 
amplitude and phase in pattern= — pi2. Similarly, for a phase modulating screen 

Amplitude in pattern =^ (specular signal) + ^2 (^j^) 
Phase in pattern=^£'i (a;,?/) /(specular signal) 
and also 

Correlation between phase in pattern and at screen = pii 
Correlation between phase and amplitude in pattern = pi2. 

It has therefore been established that the correlation between the amplitude and phase in a 
diffracted signal is zero at large distances, but at smaller distances is positive for a phase modu- 
lating screen and negative for an amplitude screen. This occurs at a distance from the screen 
given by a^lj or a^^D^/X. This is closely related to the Rayleigh distance for an optical 
system of aperture equal to the structure size of the screen irregularities. 

This result has been derived by setting up the wave solution. It may be partially con- 
firmed by a physical-optical argument, as follows. Consider a one-dimensional phase-modu- 
lating screen, with a phase profile </)(x). The signal at a point near the screen is determined 
only by the properties of the screen in its immediate neighborhood. In fact, any small portion 
of the screen centered about Xq acts as a positive lens of focal length 

The plane wave incident on the screen can be imagined as coming from a point source of 
waves, very far away. The image formed by the effective ^lens^' of the phase screen appears 
to have angular magnification=(l — s/i^)~^ where z is the distance of the point of observation 
from the screen. The image of the source appears the same brightness as the source itself, so 
the signal power arriving at z is proportional to the angular magnification. So the signal am- 
plitude is proportional to (l—z/Fy^^^; if z is assumed much less than F, approximately 



signal amplitude oc l + z/2F=l-{\2/A7r) f^l 



on substitution from (36). Any correlation between the amplitude and phase in the pattern 
must therefore arise through correlation between (p^z) and i^—dHldx^). 

In fact, it can easily be shown that if 0(x) is a Gaussianly distributed random function, 
there is a correlation between these two quantities of [p''(0)/p'' ''(0)]^/^^ where p{a) is the cor- 
relogram of ^{x). For a Gaussian correlogram, which has been assumed throughout this 
section, this reduces to l/V^; a result which is identical with that given for the same case 
in (35). It should be emphasized, however, that geometric optics can only be used near the 
screen, where the diffraction effects have not yet developed appreciably. 
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5. Effect of a Point Source of Waves 

In all the previous sections it has been assumed that there is no systematic variation of 
phase across the random screen. However, if the screen is illuminated with radiation from a 
point source a distance z' from the screen, instead of from a plane wave source, the phase of 
the wave before passing tlu'ough the screen varies with distance. The signal emerging from 
the screen is now given by 

E,{x.,,y,)=A{x,,v,:) exp {2'Ki{z'^'-^7:l-^if^'i''l\], 

If the cone angle of the diffracted waves is small, tlie important regions of the screen are for 
^0? 2/0^^' • So the approximation 

{^z''^7:lW^"''^z'^x%\1z'-Vyl\2z' (37) 

may be used. 

Exactly as in section 3, the vector signal E{x,y,z) at a distance z from the screen can be 
found and the phase components Ei{x,y,z) and E2ix,y,z) evaluated by a relationsliip similar 
to (9). In that ecpuition, the phase of the specular component is independent of x and y since 
the specular component is a plane wave. In the present cas(^, liowever, the speculai' wave 
is spherical and so the analogous ecjuation to (9) is 

E{x,y,z) = [E,(x,y,z) + iE2{x,y,z)]cxp {2Ti[{z+zr+x'+yJ^\l\} 

and for x,y<^z, 

[{z+zr+x'+yT'^^+^' + x'l2(z+z')+y'/2{z+z'). (38) 

The corresponding equation to (10) can then be written down as 

Eix,y,z) exp {-2m [{z+zr+x'+yTV^} = {W) jjjj^ ^(^o,2/o) • exp:{27r'^ [(x-x,)s 

+ {y-yo)r+{z''+xl+yiy^'+z{l-sy2-ry2)]/\}dxody4sdr 
or 

E{x,y,z,) = (\/\')ijii A{xo,yc) exp {2ir%GI\} dx^dy^dsdr 

where 

(?=(x-a:o).9+(7/-7yo)/+(^'2+^?+2/?)^/2+2(l-.^72-rV2)-[(3+^0'+^'+?/']^/^ 

By substitution from ('M) and (38), approximately 

G={x-Xo)s+{y-y.)r+z' + xl/2z'+yl/2z'+z-zs'/2-zr'/2-z-z'-xV2(z+z')-y^ 

= (x-Xo)s+{y-yo)r+xl/2z'+yl/2z'-zr'/2-zsy2-xy2{z+z')-'i//2(z+z'). 
The phase components Ei{x,y,z) and E2(x,y,z) are given by 

E,(x,y,z) = ^{E(x,y,z)} = il/2\'^) ^f(( A{xo,yo){exp (2ie/X) + exp {-2iG/\)} dx^dycdsdA 

^^^^- 1(39) 

E,{x,y,z)^j{E{x,y,z)] = (\l2\')^^^^ ^(^o,2/o){exp (2^G/X)-exp (-2iG/X) ]dx,dydsdr J 
and the spatial spectra of Ei{x,y,z) and E2{x,y,z) are given, as before, by 

Pii^,^)= j j Ei(x,y,z) . exp {—27ri(vx+xy)l\}dxdy 

rr^^ r (40) 

^2(^x)= E2(x,y,z) ' exp {-2Tri{vx+xy)/\}dxdy. 

J J — QO ^ 
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On substituting I'roni (39) into (40) and integrating with respect to s, r, x, and y, an analogous 
expression to (30) is obtained: 

Pi{v,x)^^^ cos {t{v''+X^)z{z+z')/Xz'}- jj Aixo,yn) exp {-2Ti(pXo+Xyo)iz+z')/Xz'} dx^Vo- 

(41) 

The differences between (41) and (11) can be accounted for by supposing that there are ap- 
parent values of z and X which, when substituted in (11) and in all the results derived from it, 
just allow for the fact that there is a point source of waves at a distance z[ from the screen. 

In comparing (41) and (11), the double integral will first be considered. Both v and X in the 
exponent are multiplied by the factor {z^z')lz; so the angular spectrum must be scaled down 
in width by the reciprocal of this factor, relative to its width in the plane wave case. 

The '^generalized autocorrelation function^' of Booker, Ratcliffe, and Shinn (1950) given by 
the Fourier transform of the angular power spectrum, is invariant with distance in the plane 
wave case. With point source illumination, its width increases proportionally to (l+s/^'), 
and finally increases proportionally to z. In the case of vertical incidence of waves on the 
ionosphere z=z' , and the apparent structure size of the screen irregularities is just twice the 
actual structure size. 

These results can be appreciated physically by considering the behavior of the diffraction 
pattern as X^O, and an analogy can be drawn with geometrical optics. If an amplitude screen 
(see fig. 8) is illuminated with a point source of light at a distance s' from it, the scale of the 
amplitude variations a distance z away from the screen will just be (1 + ^/^') times the scale of 
the amplitude variations at the screen. 

The other respect in which (41) and (11) differ is in the argument of the cosine. The quan- 
tity {v'^-^X^)z{z+z^)lz' replaces (j/^+x^)s. But the quantities v and x have already been found 
to be multiplied by (1 + ^/^0 ; so ^ in (30) must be replaced by 



z{z^z') 



1 



zz 



{i+zjz'y z+z' 



-'^ett) 



where Zett is the effective value of z. In the case of vertical incidence on the ionosphere, 

z=z\ and ^eff=^/2. 



SOURCE OF WAVES 

+ 




SCREEN 



Figure 8. Diagram showing the increase in structure 
size of the amplitude variations with increasing 
distance from a point source of waves. 



PLANE OF OBSERVATION 
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It is interesting to note that 

_1 J^_l 

which is exactly comparable with the object and image distance relations for a zone plate, the 
theory of which is closely related to that given above. 

6. Effect of Oblique Incidence 

If a plane wave is incident at an angle i to a diffracting screen, the equations corresponding 
to those of section 3 can be deduced by an exactly similar method. If y axis is in the plane 
of the incident wave, we may write s=sin i+Si, where (sin~^ s— sin~^5i) is the scattering angle. 
The phase of the specular signal at a distance z from the screen is now 

(27riz cos i)/\ 
and (9) must therefore be written 

Eix,y,z) = [E,{x,y) + iE2(x,y)] exp [(2Triz cos i)/X]. (42) 

The phase of tlie scattered wave at distance z is again given by 

27^^^{l + (s-2-l)-l+(r2-l)-l}-l/2/X. 

Putting s=sin i-\-Si and expanding in powers of ^i and r, this becomes 

27ri^(cos i—Si tan i—^sj sec^ i— J/'^ cos'"^ ^)/X 

and the exponent in the resulting equation for {Ei-\-jE2), analogous to (10), is 

. r 1.1^ 

27^^A {x—Xq—z tan 'i)si-{- {y—yo)r—-sl z scc^ i—-r^z cos^ i WX- 

Comparing this with the exponent in (10), it is evident that 

(i) the origin of x has been effectively displaced a distance z tan i in the plane of propa- 
gation — namely, to the ''mid-point of the path^', in propagation terminology. 

(ii) the effective distance of the observing point from the screen has been multiplied by 
factors of sec^ i and cos^ i respectively for the x and y variations in the screen. 

To make this more easily comprehensible, consider the case of oblique incidence on a 
diffracting screen with an isotropic spatial correlation function, Gaussian in form. The 
structure size of the interference pattern is found by putting 

ai=asec^ i\ 

„ . >- ii=sec^ ^ 

a2=acosHJ 

in (28). Near the screen, from (29), 

r sec^^ ^+(2/3) sec^^ i+l \^/^ 
^' \ sec24^+2secl2^+5 / 

..^ r sec^^'^+(2/3) sec^2^+l \^^^ 
' \ 5sec24i+2sec^2^+l J 

Another important effect of oblique incidence, which does not enter explicitly into this 
analysis, is the following. If the random screen is composed of a ''buckled specular reflector'^, 
it will act as a phase screen for normally incident waves. However, at oblique incidence, the 
phase change it causes in the wave is scaled by a factor cos i. This is related to the familiar 
phenomenon that an optical surface becomes more nearly a specular reflector at grazing inci- 
dence. The same effect is observed for radio waves. 
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7. Conclusion 

Detailed calculations have been given for the statistics of a radio wave diffracted by a 
random screen. The results, some of which have been used in a previous paper [Bowhill, 1957], 
now take account of anisotropy of the diffracting screen, and oblique incidence of the waves 
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